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I. INTRODUCTION 


According to the proved singularity theorems [Q the space-time singularities are the most generic 
features of Einstein’s equations. On the other hand, the nature of the space-time singularity is model 
- dependent, and still no dehnite answer about its most generic type exists. 

It is believed that the mixmaster type singularity 0, which is space-like, local and oscillatory 
can pretend to be a generic one in the general cases of a gravitational collapse and in spatially 
homogeneous cosmological models |^. The numerical works (see [Q and ref’s therein) support this 
statement for some spatially inhomogeneous cosmological models too. 

However, the recent studies of an internal structure of spherically symmetric non-Abelian 
SU (2) Einstein-Yang-Mills (EYM) black holes exhibit a new rather unexpected type of the corre¬ 
sponding singularity, which is space-like, inhnitely oscillating but not of the mixmaster type. This 
behavior of the metric is caused by the nonlinear nature of a source (Yang-Mills) held in strong-held 
region near the black hole singularity. Thus, the generic singularity inside non-Abelian EYM black 
holes can be a possible alternative to the mixmaster one if a nonlinear self-interacting matter held 
is included. 

The black hole solutions in the EYM model are very interesting objects for several reasons P|. 
They violate the naive no-hair conjecture and exhibit a discrete structure for an external solutions 
which come from the corresponding singular boundary-valued problem, imposed in a region between 
the event horizon and the spatial inhnity. Being considered dynamically, regular Bartnik-McKinnon 
solitons 0 (limited cases of EYM black holes in a limit of a shrinking event horizon) are found to 
be unstable both linearly and nonlinearly [|Tl[]. The corresponding non-Abelian EYM black hole 


solutions in an external region are also unstable under small linear perturbations, and there exist 
strong evidences that they are unstable nonlinearly [|T2 . 


The goal of the present work is to study the evolution of small but nonlinear perturbations, arising 
and propagating in the internal region of non-Abelian EYM black holes towards the singularity. We 
solve the full system of self-consistent EYM evolution equations using some kind of an adaptive mesh 
rehnement method for numerical simulations. 

The dynamics of small perturbations in black hole interior regions were investigated hrst for the 
Reissner-Nordstrom black holes ||T3| , p!^ . The qualitative predictions of an unbounded growth of 
perturbations near the Cauchy horizon were conhrmed later in a rigorous self-consistent analytical 
approach by W. Israel and E. Poisson |I3[ . 

Our investigations show that small perturbations which evolve inside non-Abelian EYM black holes 
of the generic type do not grow unboundedly and it allows us to put forward the conjecture, that 
unlike to EYM black hole solutions in an external (weak held) region, the corresponding generic 
(oscillating) internal solutions are stable, while an exceptional Schwarzschild (S) - type internal 
solutions transform to a generic one under the inhuence of nonlinear perturbations. Thus, the 
generic (oscillating) type of the space-time can be a hnal stage of a spherically symmetric collapse 
of the Yang-Mills held in an internal region of an acquired EYM black hole. 

Recently Choptuik, Chmaj and Bizoh [^] have studied the collapse of the self-gravitating YM held 
[see also 0] . They have investigated the external area of the collapsing matter up to the horizon 
formation. To penetrate under horizon it is necessary to use some kind of null coordinates in 
order to give a hnal answer on the question about the nature of the resulting space-time singularity. 

In Section II we write down the full system of EYM equations and discuss imposed initial conditions 
and background conhgurations; in Section III we briehy describe our numerical algorithm; in Section 
IV we discuss the evolution of perturbations inside generic (oscillating) EYM black hole interiors, 
and in Section V - inside S-type interiors; Section VI contains the conclusions. 
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II. THE MODEL AND FIELD EQUATIONS 


We use a spherically symmetric purely magnetic SU{2) ansatz for the Yang-Mills (YM) connection 


A = lf(t,r) - l]{T^dO - TesinOdcj)), (2.1) 

where and Tq are spherical projections of the SU{2) generators. The general ansatz for the 
spherically symmetric SU{2) YM held admits also the second independent function, originated from 
the Aq component of connection, but we omit it here. 

The four-dimensional spherically symmetric metric tensor also admits two independent functions. 
Therefore we can choose the following parametrization of the interval: 

ds^ = dt^ — ^dr'^ — r‘^d9‘^ — sin^ OdcfP'. (2.2) 

A 

Both metric functions A and a as well as the YM function / depend on r and t variables. 

Let us denote 


D = —, n = -^, (j) = f, {f = dr, ' = dt). 

r Da 


After that the full set of EYM equations looks as follows: 


1 — D' = K —-|- K — 
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- = -(02 
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D = -2k 
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aU(j)D^ 
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r = 0, 
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\ r 
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n = a 

■0 = 0 . 


f{p - 1) 


(2.3) 

(2.4) 

(2.6) 

( 2 . 6 ) 

(2.7) 

(2.8) 

(2.9) 


Note, that equation 


corresponds to the tt component of the Einstein equations Gl = nTj:, 


equation ( ^.41) - to the difference of tt and rr components Gl — = k{T^ — equation (^.5[ ) is 

the rt component G^ = kT[ of the Einstein equations, equation (|2.6|) is the dehnition of fl, equation 
(EZP is the dehnition of 0, equation (|2.^ ) is the Yang-Mills equation of motion, and equation (13 
is the requirement of / to be smooth: ~ 

It is well-known that the Schwarzschildean radial coordinate r becomes the temporal coordinate, 
while t becomes the spatial one in the region under the event horizon of a black hole and the 
dynamics in the black hole interior region is described by the evolution equations along r, together 
with constraints at each slice r = const. 

Now (p7^) , ( p.4|) and (p^) are evolution equations, (p78|) and ( p.9|) are wave equations, whereas (p.5|) 
and (|2.6| ) occur to be conserved constraints, which hold automatically. Indeed, after diherentiating 


( p.6|) with respect to r one has identically 


'Da. 


n)Y|./'-(—n| =^,-,# = 0. 

r / at \ r 
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Following the same lines, if we denote the relevant combination in (|2.5|) as 


'-) = i) + 2 k 






then it can be easily shown that 

7 =7—• 

cr 

So, if 7 = 0 at the initial surface r = const then 7 will be zero during the evolution along r. As a 
result, there are no dynamical constraints in our EYM system in the black hole interior region; both 
constraints are kept during the evolution automatically and the system is effectively described by 
the equations (^), (13), (g3), (13) and O. 

However, to realize the numerical algorithm, a little different representation of unknown functions 
turns out to be more effective. Let us introduce an auxiliary field 


X = Aa/r (2.10) 

as a dynamical variable instead of A. Unlike A, the held y does not exhibit oscillations for the 
background generic solutions in the interior regions of EYM black holes; in terms of y the local 
speed of light (the slope of the characteristics) is equal to Y 
Then, approaching the black hole singularity at r —> 0, it is more suitable to use an inverse 
coordinate R= 1/r for the numerical integration. 

Now the resulting set of unknown functions, used for numerical study of our PDE system consists 
of 


X = AaR, a, ct> = xR^ n = / 


(here and below f = Or ), and the full system of equations obtained from (^.3|), (^.41), (p^.S]), (^.9]) 
hnally has the form: 


X' ^ . (f (7 -1)= - j ,). 


a 

a 


K 


R^X- 


:(0^ + n2), 


^ ^ ^ ^ = afif - 1 ), 


dt \R^x, 


dt \R^x 


/' = 


R^X 


Apart from these equations we also use the constraints 

n = /, 
l = + 


X 


X 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.17) 


in order to set the initial data and to keep the control of the accuracy. 

We set the coupling constant k = 1 hereafter without loss of generality. 

After the completion of numerical calculations we display the results again in terms of metric 
function A(r, f) to get more insights on their physical meaning and to compare them with the 
corresponding background configurations. 
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A. Background configurations 


We study the evolution of perturbations, propagating on a homogeneous (t -independent) back¬ 
ground configurations which correspond to the interiors of spherically symmetric EYM black holes. 
These configurations are the solutions of the system 


V r / \ r 


= 1 


iP - 1 )^ 


(2.18) 


(In a)' = 


2 / 


l2 


obtained from (|2.8|) , (|2.3|) and (|2.4|) , neglecting the t -dependence in the domain r < Vh- The initial 
conditions are imposed in the small vicinity of the surface r = which is implied to be the position 
of a simple (not double) event horizon. As it was shown recently the generic solution of the 

system ( p.l^ ), corresponding to the interior of a static spherically symmetric EYM black hole, has no 
Cauchy horizons, and the metric exhibits an infinitely oscillating behavior (but not of the mixmaster 
type) with an amplitude, unboundedly growing towards the space-like singularity. 

The oscillating structure of the metric for the generic solution originates from the features of 
the corresponding 2-dimensional autonomous dynamical system. This system effectively describes 
generic solutions in the regime when some irrelevant terms in (|2.18|) (1 is negligible in comparison 
to {p — l)^/r^ and the YM function / is set equal to constant, / = const, while drf ^ 0 remains 
dynamical ) are neglected near r —> 0 P|. In the interior of the EYM black hole the metric passes 
through an infinite series of “almost Cauchy horizons” in the maxima of oscillations which alternate 
by subsequent huge falls of the metric function A in the minima; the frequency of oscillations of the 
metric exponentially grows as the singularity is approached. The approximate recurrence formulas, 
obtained in §], allow one to describe the behavior of the EYM system in such a “strong oscillation” 
regime with an accuracy, improving towards the singularity. 

However, for a typical generic EYM internal black hole solution the “strong oscillation” regime 
described above does not start just in a vicinity of the event horizon. Depending on the initial 
conditions on the event horizon, before the first huge fall of the “strong oscillation” regime begins at 
some r = Tosc.^ the solution is determined by the complete system (|2.18|) with all terms relevant in 
the significant domain rose. < r Kr^. In this domain the metric function A, being negatively defined, 
also can admit a few oscillations with a small amplitude (in comparison with “strong oscillations”); 
we call this regime a “weak oscillation” regime. 

In the present paper we consider the evolution of perturbations starting from an internal vicinity of 
the event horizon and then propagating through the “weak oscillation” region and the first huge fall 
of the metric function A in the “strong oscillation” regime up to the next “almost Cauchy horizon”. 

As it was also shown in 0.0. for some discrete values of initial parameters, EYM spherically 
symmetric black holes also admit the standard Schwarzschild and Reissner-Nordstrom (RN) interiors 
(third, so called HMI internal solution is not asymptotically fiat and we do not consider it here); 
however such interiors are rather exceptional cases and they form only a subset of a zero measure 
in the space of all EYM internal black hole solutions. The evolution of perturbations inside S -type 
EYM black holes is investigated as well and the results will be discussed below. 
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B. Initial conditions 


The Cauchy problem for the system of equations ( p.ll|) - (|2.15|) can be set as follows: we set 
/, 0, X at some space-like surface r = const and define If and a from (|2.16|) , (|2.17|) . Once 


being imposed, the constraints (|2.16|) and (p.l7|) will be satisfied during the evolution along r due to 
the equations ( 2.11 ) - (|2.15|) . The latter can be solved using a standard finite - difference numerical 
technique. 

Note, that we suggest the initial space-like hypersurface r = const < Vh situated with the necessity 
under the external event horizon. Therefore the perturbations given at this hypersurface have a 
status of fluctuations, acquiring in an internal region of a black hole and they are not related 
anyhow with those, propagating inward from an external region through an event horizon. 

The initial values of / and 0 can be set independently. Straumann and Zhou [0 considered two 
types of initial perturbations for the case of Bartnik-McKinnon regular solitons and the external 
region of EYM black holes. They called the perturbations class I if only the YM function / is 
initially perturbed, while the time derivative of the perturbation equals to zero 50 = 0; it means 
that the initial kinetic energy of the YM held vanishes. In class II perturbations the function / 
remains unperturbed, but 50 and thus the kinetic energy does not vanish initially. 

We used both these classes of initial perturbations for our problem in the internal region of EYM 
black holes. The considered initial data come from the asympthotics near the external event horizon. 
Since the horizon is not a regular point, we cannot set all Cauchy data independently. The data 
must satisfy the series expansions originated from the Einstein-Yang-Mills equations. 

The hrst order of asympthotics for background t -independent EYM equations give A = di{r — 
Th), a = ao + CTii^r — Th), f = fo + — Vh), where Vh is the radius of the event horizon and /o, 

(To are free parameters. Other parameters are expressed in terms of these free parameters as follows: 


<ii = >■» - (/o - /i =/o(/o - CTi = 2(T|,/,7ri. 


(2.19) 


Thus, in order to define the background solution, we should set the values /o, (Tq ; the coeffi¬ 
cients di, fi and ai are determined from (|2.19| ). 

We choose an initial surface at the distance h (one r - step) under the horizon. The equilibrium 
(background) values of / and y at this surface are / = /o — /ih, y = —dia^h/rh- 
Class I perturbations. In the case when small, for example. Gauss - like perturbation 




( 2 . 20 ) 


is added to the equilibrium YM function /, then the held 11 is determined by the equation ( p.l6| ). 


i.e.. 


Il{t)\r=rf,-h = -2Ks{t - to)e 




For initial y we use the equilibrium value and initial 0 is also set to be unperturbed. Thus, class I 
initial perturbations are dehned as follows: 


X{t)\r=ru-h = -diaoh/rh, 
(p{t)\r=r^_h = -diaofih/rl 


( 2 . 21 ) 


with Yl{t)\r=rh-h determined above. 

Class II perturbations. In this case the YM function / remains initially unperturbed along with 
the corresponding function 11, while the function 0 gets independent deviation from the background: 
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f{t)\r=r,,-h = fo-flh, 

U{t)\r=r^_h = 0 , 

Xit)\r=r^-h = -diaoh/rh, 

(j){t)\r=rf^_h = -diaofih/rl + Ke-’^d-to? _ 


Both classes of initial perturbations give rise to two scattering waves. One of them propagates 
with time r in a positive spatial t -direction (“ingoing” wave) and another one propagates in the 
negative t -direction (“outgoing” wave). After some interference in a region around t = Iq they move 
independently from each other in opposite spatial directions. We investigated the evolution of class 
I and class II initial perturbations on various EYM black hole interior backgrounds and found their 
behavior to be the same after the “ingoing” and “outgoing” waves are completely divided in the 
space. Moreover, if we set the initial perturbation with a symmetry with respect to the center of 
perturbation to, this symmetry is conserved during the evolution. 

That is why one can consider without loss of generality the third class of initial perturbations 
in addition to the described above; the initial perturbations of this class determine only one initial 
wave, propagating, say, in a positive t -direction. This class occurs to be a combination of class I 
and class II and is dehned by the following way. 

One-wave class of initial perturbations. Perturbations of this class are the same as in class I for 
the functions / and 11. In addition, the perturbation of (j) is chosen not independently, but in such 
a way, that the initial “outgoing” wave is canceled: 

fit)\r=r,-h = + 

Il{t)\r=r^-h = -2Ks{t-to)e-^d-to)\ 

x{t)\r=rH-h = -diaoh/vh, ^ ^ 

(j){t)\r=r^_h = -diaofih/rl-2Ks{t-to)e-^^^-^°'>\ 


We will use this class of initial perturbations for the illustration of their nonlinear evolution in 
subsequent Sections. 

Now the initial value for the held a for all considered classes of perturbations can be determined 
by the numerical integration of the equation (|2.17|) . We integrated it from the left to the right, since 
the unperturbed functions are situated at the left side at f << fo¬ 


ul. NUMERICAL METHOD 


The specihc feature of the strong - held dynamics requires a numerical code which must be able 
to resolve a solution on extremely small scales. Usually some kind of an adaptive mesh rehnement 
algorithm is implemented for these problems. 

We controlled the resetting of grids “manually”. The accuracy was determined by means of 
calculation of the constraint When the absolute value of this constraint at least at one point 

became large than 10“^ the program saved its state and stopped. Then we investigated output data 
visually, determined the “bad” region and set new grid steps. Starting again, the program picked 
out necessary region and recalculated the data according to new grid steps. In order to set 11 and a 
according to new data, the constraints (|2.16|) , (|2.17|) were also recalculated. 


A. Algorithm 


The system of equations (p.ll|) - ( p.lSj) was solved by means of MacCormac predictor - corrector 


scheme (see, for example, J^). We tested many modihcations of this scheme and convinced ourselves 
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that in this case various hnite differences used for forward/backward difference lead approximately 
to the same accuracy. The results can vary from one domain to another, but not very strongly. The 
most universal approach appears to be the conventional MacCormac forward/backward difference, 
so we used it for our calculations. 

The constraint ( p.l6| ) does not contain derivatives of unknown function, so we integrated it using 
a simple midpoint formula: 


fi —1 

ti—i 

For integration of the constraint ( p.l7|) we tried to use several methods, such as various Runge- 
Kutta schemes, predictor-corrector methods, etc. The best result was obtained by the simplest 
midpoint formula: 


O'i+l — CTi-l + {ti+1 ~ ti_i) 


(Ti 


Xi 


f Xi+1 - Xi-1 

\ ^i+1 ^i—1 




(3.1) 


This result is obvious because we used centered differencing for monitoring of the constraint ( p.l7| ). 

We calculated the maximal slope of the characteristics on each slice R = 1/r = const and 
then chose h/i (“time” R - step) according to the Courant criterion, hji = htC/2. We denote the 
“spatial” current integration step by ht and the maximum of the slope of characteristics on a current 
slice as C = max(yi?^). 

The value of C tends to cxo as r —> 0 (or R oo). In order to prevent a destruction of the 
numerical algorithm, we bounded it by some value Cmax- Surely, the Courant criterion was satisfied 
with guarantee for all Cmax < C. The final expression for looks as follows: 


hn — 


Cht 
2 ’ 
Cmax ht 
2 ’ 


if C < Cmax 
if C > Cmax 


The value of Cmax was set “by hand” and varied from 10 to 10^ in different regions. 


B. Realization remarks 

The algorithm was realized as ANSI C program. For real numbers the type double was used. 

The initial “spatial” (f) step was set equal to 10“^, the relevant initial interval was —30.0 < t < 
30.0. So, the number of one slice points we started with was about 6000. During the calculation, 
this number may increase up to 10®. In output the program saved only the part of the points, and 
the distribution function of the saved points was defined according to the estimation of an error, 
produced by the usage of the coarse grid. 

The initial data were calculated in the coordinate r and then converted into the coordinate R. 

The typical time needed for evaluation of one configuration was about 50 hours in Pentium-166. 

IV. GENERIC SOLUTIONS 

We have tested various background generic internal configurations with values of Vh in a range 
from 10“^ up to 10®, corresponding to both asymptotically flat (black hole) and not asymptotically 
flat EYM solutions in an external region. Similarly to the background configurations, the evolution 
of perturbations under the event horizon turns out not to be sensitive to the background type in 
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the external region (either asymptotically flat or not) and looks qualitatively the same for both 
considered types. 



-0.688135 -ln(r) 12.951982 

FIG. 1. Background A under the event horizon. One can see two oscillations of the function A in 
the “weak oscillation” region (second minimum at — Inr = 3.81671 is non-observable in this scale) 
and the first oscillation in the “strong oscillation” region. Although the first fall of ^background is 
not very huge, it is described, however, by the reduced system of equations with a good accuracy 
and therefore corresponds to the “strong oscillation” regime. Thick grey lines below mark regions 
depicted in 3-D plots of 5A: a) corresponds to Fig.|, this region is very small; b) corresponds to 
Fig.H; c) corresponds to Fig.^ and d) corresponds to Fig.^. 


We were able to follow the evolution of perturbations of both classes (I and II) with different (small 
enough) initial amplitudes, propagating through the “weak oscillation” region and the first huge fall 
of the function A up to a vicinity of the first “almost Cauchy horizon” of the “strong oscillation” 
region. The common result is that during propagation towards the singularity the considered per¬ 
turbations do not grow up unboundedly and look like a light ripple on the background solutions. 
Therefore the strong numerical evidence in favor of nonlinear stability of generic (oscillation type) 
internal EYM black hole solutions is obtained. 

To illustrate this statement we have chosen an appropriate background configuration. Being 
asymptotically not flat in exterior, it exhibits all important features of EYM black hole interiors 
and allows one to plot both “weak” and “strong oscillation” regions in the same figure with enough 
resolution. The chosen background is characterized by the following parameters 

rh = 2.{) (72 = 0.5, - Inr =-0.69314718); 

/o = -0.302072; 
cTo = 0.9, 

and the corresponding curve of the metric function A is plotted in Fig. |l[ 

For the illustration purposes we have used the one-wave class of initial perturbations. According 
to Section II, it is determined by a f -dependent deviation of the YM function / from its background 
value which induces the perturbation of 11; the perturbation of the function 0 is defined to cancel 
the initial wave, propagating to the negative t -direction. 

We plot below the evolution of perturbations with initial parameters in (|2.23| ): 

s = 100.0; to = -10.0 

for two different amplitudes Ki = lO”"^ and K 2 = 10“^ in (|2.20|) (we will call them “small” and “big” 
perturbations respectively in the further discussion). These perturbations induce deviations of all 
other functions according to the full set of EYM equations and the resulting perturbation looks like 
a nonlinear wave, propagating in the positive spatial t -direction with time r, directed to r —>■ 0. 

To show the relative size and shape of perturbations one can normalize the corresponding functions 
by their equilibrium values. The value SA 2 , obtained as 6 A 2 = A 2 — Abackground for “big” initial 
perturbation, is plotted in four figures: Fig.||, Fig.|], Fig.| and Fig.|[ The considered r - regions are 
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shown in Fig.|l|by thick grey lines. The plots illnstrate a strongly nonlinear natnre of the evolntion: 
the pertnrbation of the metric fnnction A changes the relative sign and the shape as singnlarity 
r —> 0 is approached. The pertnrbations of other relevant fnnctions behave in a similar way and we 
do not plot them here. 


-ln(r) 



-12 

FIG. 2. The beginning of the evolntion of pertnrbation of 6 A 2 = A 2 — Abackground- The initial 
pertnrbation looks like a step. The folds of this step remain to be of a constant height dnring the 
beginning of evolntion. Further the left fold (greater t) will change while the right is constant by 
definition. This figure corresponds to region a) in Fig.|l|. The arrow points to the direction of the 
evolution. 



.5 

. 0015 


t 

FIG. 3. Plot of 6 A 2 in the beginning of the “weak oscillation” region. This figure corresponds to 
b) region in Fig.|^ The arrow points to the direction of the evolution. 
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FIG. 4. Plot of 6 A 2 in the “weak oscillation” region. This figure corresponds to c) region in Fig.|^. 
The arrow points to the direction of the evolution. 



FIG. 5. Plot of 6 A 2 in the “strong oscillation” region. This hgure corresponds to d) region in Fig.|l|. 
The arrow points to the direction of the evolution. 

However, the perturbations do not grow unboundedly and therefore they do not destroy the back¬ 
ground solution. It is convenient to plot the absolute values of the maximal deviation of all relevant 
functions from their backgrounds as perturbations are developed. The corresponding plots are rep¬ 
resented in Figures 1-0 for both “small” {K = 10 and “big” {K = 10 initial perturbations 
along with the background functions. 

These plots exhibit most general features of the perturbations, propagating in the different con¬ 
sidered generic backgrounds: the evolution of “small” and “big” perturbations looks very similar for 
all relevant functions, so the difference in initial amplitude does not produce differences in shapes 
and characters of perturbations; the amplitude of perturbations of A is approximately proportional 
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-0.688135 -ln(r) 13.137556 

FIG. 6. A background (bottom), absolute value of the maximal deviation from the background 
(middle - “small” initial perturbation, top - “big” initial perturbation). 



FIG. 7. A in the “weak” oscillation region: background (bottom), absolute value of the maximal 
deviation from the background (middle - “small” initial perturbation, top - “big” initial perturba¬ 
tion). 
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FIG. 8. a: background (bottom), absolute value of the maximal deviation from the background 
(middle - “small” initial perturbation, top - “big” initial perturbation). 



-0.688135 -ln(r) 13.137556 

FIG. 9. /: background (bottom), absolute value of the maximal deviation from the background 
(middle - “small” initial perturbation, top - “big” initial perturbation). 

to the absolute value of A itself (it grows in the “strong oscillation” regime following the fall of 
^background and then decreases as ^background approaches the “almost Cauchy horizon”); perturbations 
of all other functions demonstrate bounded nonlinear oscillations in the “weak oscillation” region, 
then they become almost constant ones in the “strong oscillation” region up to a vicinity of an 
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’almost Cauchy horizon”. 



-0.688135 -ln(r) 13.137556 

FIG. 10. 0: background (bottom), absolute value of the maximal deviation from the background 
(middle - “small” initial perturbation, top - “big” initial perturbation). 



-0.688135 -ln(r) 13.137556 

FIG. 11. fl; absolute value of the maximal deviation from the background (bottom - “small” initial 
perturbation, top - “big” initial perturbation); the background value of If is zero by dehnition. 

We have followed the evolution of the perturbations towards the r = 0 singularity up to the first 
“almost Cauchy horizon” in the “strong oscillation” regime. However, there are no physical reasons 
to consider the next huge oscillation of the metric function A, since in this region the magnitude of 
the Riemannian squared scalar exceeds the Planckian value > ^/L^i) and the classical 

description of space-time is no longer valid overthere. So, the only problem remains to penetrate 
numerically through the first “almost Cauchy horizon” in the “strong oscillation” regime. We are 
going to consider this problem separately along with RN -type using the more precise numerical 
code. 
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V. SCHWARZSCHILD - TYPE SOLUTIONS 


Schwarzschild and Reissner-Nordstrom types of internal backgronnd EYM black hole solntions are 
of exclnsively types since they form a set of zero measnre in the space of all initial data , 0] . They 
exist only for some discrete valnes of {vh, /o) on a horizon and any internal backgronnd solntion of 
( P.18| ) with an arbitrary small deviated initial data from these discrete valnes does correspond to a 
generic (oscillating) type. 



t=20 t=10 t=0 t=-10 t=-20 



FIG. 12 . The wide table-like pertnrbation imposed at the initial hypersurface came from the 
horizon at t = -|-cx) moves in the negative t -direction (grey region). In the r, t coordinates a) the 
perturbation almost stops its motion in t -direction and only goes with time r to the singularity 
r = 0. It will never cross t = —20 hypersurface, but it crosses an arbitrary section t = const > 0 
in the past with necessity, as it follows from the conformal diagram, depicted in b). The top of 
the “table”, which corresponds to the cross - section t = 0 can reach the singularity r = 0 before 
perturbations from the edges of the initial “table” (curve A) will reach it. 

It can be expected that small t -dependent perturbations added to the initial data of these excep¬ 
tional interiors will produce the transformation to some generic type during the nonlinear evolution 
towards r —> 0. We have investigated the dynamics of perturbations on various Schwarzschild - type 
EYM black hole internal solutions and convinced ourselves that this transformation really takes 
place and therefore S-type interiors are occurred to be unstable. 

Indeed, according to (|2.21|) , ( p.22| ) and ( p.23| ), the perturbations produce deviations from the 
background value of the YM function / in the leading order of series expansions with respect to r 
near the event horizon; it is true by dehnition for class I and one-wave class of initial perturbations, 
while the initial perturbations of class II produce similar deviation of the YM function / at the next 
step on r as perturbation starts to evolve. 

Thus, only if waves, produced by the initial perturbation, are suppressed fast enough during their 
evolution with time r directed towards r —0, there is a chance for the exceptional internal solution 
to conserve its S-type. However, our numerical studies show that perturbations are not suppressed 
during their evolution inward S-type EYM black hole interior and, as a result, S-type singularity 
transforms to the generic oscillating type. 

To illustrate this process it is convenient to investigate S-type interior, perturbed by the initial 
table-like “outgoing” deviation from the background, which produces the corresponding wave, prop¬ 
agating with time r from t = -|-oo to f < 0. We integrated the equation (p.l7|) from the left to 
the right, the region with large enough negative t corresponds to the originally unperturbed S-type 
solution since the perturbation never can reach this region [see Fig-Hl. 

Investigating the small r -vicinity of the singularity in different spatial t -regions one can answer 
the questions about the type of the resulting solution. 

We chose the shape of the initial perturbations of / to be a Gauss - like curve, broken at the top 
by large table - like insertion [see Fig.p!^: 
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FIG. 13. Initial / perturbation for the Schwarzschild - type solution. 


r t < toi, 

5f{t) = \ K, toi <t < to2, 

[ t > to2. 

Values of the parameters for the initial data were set equal to: 

Th = 0.613861419, 

/o = -0.8478649145, 
do = 0.289427236, 

K = 0.02, 
foi = —10.0, 
to2 = 10.0. 


This internal background S-type conhguration corresponds to the asymptotically flat (black hole) 
solution with iV = 1 (nodes of YM function) in the external region. We choose the value of 0 in the 
way that the outgoing wave goes from t = +cx5 to f < 0: 

( -diaofih/rl + 2Ks{t - , t < toi, 

0(t) = } -diaofih/rl, toi<t < to2, (5.1) 

[ -diaofih/rl + 2Ks{t - t > to 2 - 


This shape of the initial perturbations permits us to investigate three different spatial cross - 
sections t = —20, t = 0 and t = 20 independently, since the perturbations from the edges of the 
initial “table” (at t = —10 and t = 10) can not reach these points during the evolution [see Fig.lK . 

The solution in the hrst spatial cross - section corresponds to the background Schwarzschild - type 
solution. Initial perturbations cannot reach this region, and we obtain the typical Schwarzschild - 
like behavior [see Fig.|T^ a) and b)]: the function y goes to the constant, approaching the singularity 
r = 0 [Fig.|T^ b)]. 

The generic - like y goes to 0 as r —>• 0. It is easy to see that the cross - section through the top 
of the “table” at t = 0 demonstrates just this type of the behavior [FigJT^ c)], since the considered 
perturbation are not suppressed during the evolution towards r —> 0 and the singularity becomes of 
the generic (oscillating) type. 
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FIG. 14. Various cross - sections of x for the Schwarzschild - type internal solution. The cross 
- section t = —20 (a) corresponds to the non-perturbed background solution: approaching the 
singularity r = 0, the function y goes to the constant (b), while for any generic (oscillating) internal 
solution the function y goes to 0 as r —> 0. The cross - sections t = 0 and t = +20 (c), (d) 
demonstrates that the internal solution becomes the generic one. 

Considering a small vicinity of r = 0, one can see that the resulting solution has already the generic 
type [Fig.|^ d)] also in the third spatial cross - section t = +20. The transformation of S-type to 
the generic one in this region is caused by the shift of the apparent horizon position. The pulse of 
the considered perturbation shifts the ADM mass of the system and the position of the apparent 
horizon ^ Vh + while the relevant value of the YM function on the apparent horizon /q 
remains unperturbed. As a result, new initial data {rh + Svh, fo) for the homogeneous system (p.l8|) 
now corresponds to the generic type of the internal solution. 

VI. CONCLUSIONS 

We have investigated the dynamical evolution of small initial perturbations in space-time re¬ 
gions which correspond to internal part of a spherically symmetric black hole in non-Abelian purely 
magnetic SU{2) EYM theory. The obtained results give a strong numerical evidence in favor of 
nonlinear stability of the generic (oscillating) type of EYM black hole interiors while the excep¬ 
tional Schwarzschild - type interiors turn out to be unstable and transform to the generic type as 
perturbations are developed towards a singularity. 

Now one can expect that the generic (oscillating) type of the EYM black hole singularity is stable 
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as well with respect to the perturbations, penetrating into the internal region from the exterior 
through the event horizon. Moreover, the generic interior solution can pretend to be the hnal stage 
of a spherically symmetric collapse of the Yang-Mills held after the event horizon is formed. To 
check this expectations one should use the null coordinates and the more precise software tools to 
attack the problem numerically. This work is in progress now. 
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